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ABSTRACT. We consider the nonlinear Schrodinger equation with cubic (focus- 
ing or defocusing) nonlinearity on the multidimensional torus. For special small 
initial data containing only five modes, we exhibit a countable set of time layers 
in which arbitrarily large modes are created. The proof relies on a reduction to 
multiphase weakly nonlinear geometric optics, and on the study of a particular 
two-dimensional discrete dynamical system. 



1. Introduction and main result 
We consider the nonlinear Schrodinger equation 

(1.1) id t u + Au = X\u\ 2 u, x G T d , 

with d ^ 2, where the sign of A G {— 1, +1} turns out to be irrelevant in the 
analysis below. In the present analysis, we are interested in the description of some 
energy exchanges between low and high frequencies for particular solutions of this 
equation. We will consider solutions with small initial values: 

(1.2) u(0, x) = 6uo(x), 

where u G ff 1 (T d ) and < S <C 1. Replacing u with 5~ l u, (TTTTb— (fTT2l> is 
equivalent to 

(1.3) id t u + Au = Xe\u\ 2 u ; u(0,x) = uq(x), 

where e = 5 2 . Viewed as an infinite dimensional dynamical system in terms of 
the Fourier variables of the solution, such an equation is resonant in the sense 
that all the eigenvalues of the Laplace operator are integers only, making possible 
nontrivial vanishing linear combinations between the frequencies of the linear un- 
perturbed equation (e = 0). In such a situation, the perturbation theory cannot be 
directly applied as in |[T1 [2 [3] [9] [101 LT21 Ll5]| . Let us recall that in all these works, the 
Laplace operator is perturbed by a typical potential making resonances genetically 
disappear. In such situations and when uq{x) is smooth enough, it is possible to 
prove the quasi preservation of the Sobolev norms of the solution over very long 
time: polynomial (of order e~ r for all r) as in 0, exponentially large as in |[T3l . 
or arbitrary large for a set of specific solutions as in ifTOll . 

In the resonant case considered in this paper, there is a priori no reason to ob- 
serve long times bounds for the Sobolev norms of the solution. Despite this fact, 
Eqn. (11.11 ) possesses many quasi -periodic solutions (see B1I201). 
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On the other hand, it has been recently shown in that in the defocusing case 
(Eqn. (11.11) with A = 1), solutions exist exhibiting energy transfers between low 
and high modes which in turn induce a growth in the Sobolev norm H s with s > 1. 
Strikingly, such phenomenon arises despite the fact that L 2 and H 1 norms of the 
solution are bounded for all time. 

The goal of the present work is to describe quantitatively such energy exchanges 
in the case of a particular explicit initial value uq{x) made of five low modes. Since 
we work on T d , the solution u takes the form 

where Uj(t) G C are the Fourier coefficients of the solution, and as long as t 
does not exceed the lifespan of u. Here, for j = (ji, ■ ■ ■ ,jd) £ ^ d an d x = 

(xi, . . . ,x d ), we have j-x = jixH VjdXd- We also set \j\ 2 = j 2 ^ hjj G N. 

Let us introduce the Wiener algebra W made of functions / on T d of the form 

/(x) = y, '>,< 

such that (bj) jeZd G £ 1 (Z d ). With this space is associated the norm 

ll/lk = EN- 

Our main result is the following: 
Theorem 1.1. Let d ^ 2, and u G C°°(T d ) given by 

Uq{x) = 1 + 2COSX1 + 2COSX2- 

For A G {±1}, the following holds. There exist sq, T > such that for < e ^ £q, 
(11.31 ) has a unique solution u G C([0, T/e]; W), and there exist Cq > and C > 
such that: 

Vj G AC, 3 Cj - 7^ 0, Vt G [0,T/e], Uj- (t) - c^rt) 1 ^ 2 " 1 < (C et) li|2 + Ce, 
where the set AC w given fry 

(1.4) A/"* = {(0, ±2P), (±2P, 0), (±2P, ±2P), (T2 p , ±2?), p G N} x {0 zd - 2 }. 
Arbitrarily high modes appear with equal intensity along a cascade of time layers: 

V 7 G]0,1[, V0< -, Va>0, 3eiG]0,e ], Vee]0,ei], 



Vj G AC, |j| < a ( log- 



2 



i-7/(lil 2 -i) 



> — . 
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This result expresses the possibility of nonlinear exchanges in (11.31 ): while the 
high modes in the set AC are equal to zero at time t = 0, they are significantly 
large in a time that depends on the mode. As this time increases with the size 
of the mode, this is an energy cascade in the sense of Q. To our knowledge, this 
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result is the first one where such a dynamics is described so precisely as to quantify 
the time of ignition of different modes. 

The proof of this theorem relies on the following ingredients: 

• An approximation result showing that the analysis of the dynamics of (11.11) 
over a time of order 0(1/ e) can be reduced to the study of an infinite 
dimensional system for the amplitudes of the Fourier coefficients uj (in 
a geometric optics framework). Let us mention that this reduced system 
exactly corresponds to the resonant normal form system obtained after a 
first order Birkhoff reduction (see ifTTl for the one dimensional case). We 
detail this connection between geometric optics and normal forms in the 
second section. 

• A careful study of the dynamics of the reduced system. Here we use the 
particular structure of the initial value which consists of five modes gener- 
ating infinitely many new frequencies through the resonances interactions 
in the reduced system. A Taylor expansion (in the spirit of jO) then shows 
how all the frequencies should be a priori turned on in finite time. The 
particular geometry of the energy repartition between the frequencies then 
makes possible to estimate precisely the evolution of the particular points 
of the set A/"* in (11.41) and to quantify the energy exchanges between them. 

The construction above is very different from the one in (H. Let us mention that 
it is also valid only up to a time of order 0(1/ s) (which explains the absence of 
difference between the focusing and defocusing cases). After this time, the nature 
of the dynamics should change completely as all the frequencies of the solution 
would be significantly present in the system, and the nature of the nonlinearity 
should become relevant. 

2. An approximation result in geometric optics 

For a given element (aj)j € %d G £ 1 (Z d ), we define the following infinite dimen- 
sional resonant system 

(2.1) idj = A a kd£a m ; a,j(0) = atj. 

(k,i,m)elj 

where Ij is the set of resonant indices (see IfTTl ) associated with j defined by: 

(2.2) Ij = {(k,£,m) G l? d | j = k-e + m, and \j\ 2 = \k\ 2 - \£\ 2 + \m\ 2 } . 
With these notations, we have the following result: 

Proposition 2.1. Let uq(x) G W and (aj)j & z d e ^ 1 (Z d ) its Fourier coefficients. 
There exists T > and a unique analytic solution (cij)j^d : [0, T] — > £ 1 (Z rf ) to 
the system (12- 1 b - Moreover, there exists Eq (T) > such that for < e ^ Eq (T), the 
exact solution to (11.31 ) satisfies u G C([0, — ]; W) and there exists C independent 
of e G]0,£o(T)] such that 

sup \\u — v £ (t)\\w ^ Ce. 
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where 
(2.3) 



v £ (t,x) = a,j(et)e 



ij-x-it\j\ Q 



Remark 2.2. Even though it is not emphasized in the notation, the function u ob- 
viously depends on e, which is present in (1 1 -3b - 

We give below a (complete but short) proof of this result using geometric optics. 
Let us mention however that we can also prove this proposition using a Birkhoff 
transformation of (11.31 ) in resonant normal form as in ifTTl . We give some details 
below. There is also an obvious connection with the modulated Fourier expansion 
framework developed in lPT4l in the non-resonant case (see also |[T9ll ). 

2.1. Solution of the resonant system. The first part of Proposition 12. II is a con- 
sequence of the following result: 

Lemma 2.3. Let a = {oij)j eZ d £ £ 1 (1i d ). There exists T > and a unique 
analytic solution (aj)iez d ■ —> I 1 {IS 1 ) to the system (12- 1 b - Moreover, there 

exists constants M and R such that for all n £ N and all s ^ T, 



(2.4) 



Vj £ 



d n cij 

; (S) 

dt n K 1 



< MR n n\ 



Proof. In J6l, the existence of a time T\ and continuity in time of the solution 
a{t) = (a,j(t))j £Z d in t is proved. As £ is an algebra, a bootstrap argument 
shows that o(t) £ C°° ([0,T 1 ];£ 1 (Z' i )). From (O) we immediately obtain for 
se [0,7i], 

||a(s)||,i <3||o(a)||^, 

and by induction 



,(«)/ 



2n+l 



( S ) ^ <3-5---- (2n + l)||a( S )||£ 
where a( n ) (i) denote the n-th derivative of a(t) with respect to time. This implies 



(0) 



3 • 5' 



• (2n + l)||a|| £ 2 r +1 < ||a||^n! (3||a||| 1 ) n , 

which shows the analyticity of a for t < T 2 = ^\\a\\J^. The estimate (12.41 ) is then 
a standard consequence of Cauchy estimates applied to the complex power series 
En G N ^a (n) (0)-? n defined in the ball B(0, 2T) where 2T = mm(T 1 ,T 2 ). □ 

2.2. Geometric optics. Let us introduce the scaling 

(2.5) t = st, x = ex, u(t, x) = u e (et, ex) . 
Then dl.31 ) is equivalent to: 

(2.6) ied t u £ + e 2 Au £ = Ae|u £ |V ; u £ (0,x) 



u 



ie2 



In the limit e — > 0, multiphase geometric optics provides an approximate solution 
for ( 12.61 ). The presence of the factor e in front of the nonlinearity has two conse- 
quences: in the asymptotic regime e — > 0, the eikonal equation is the same as in 
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the linear case A = 0, but the transport equation describing the evolution of the 
amplitude is nonlinear. This explains why this framework is referred to as weakly 
nonlinear geometric optics. Note that simplifying by s in the Schrodinger equation 
(I2.6I ). we can view the limit e — > as a small dispersion limit, as in e.g. |[T6l . 

We sketch the approach described more precisely in |6j. The approximate solu- 
tion provided by geometric optics has the form 

(2.7) v £ (t,x)= ^ a i (t)e^( t ' x )/ £ , 

where we demand u e = v e at time t = 0, that is 

aj(0)=aj ; <^-(0, x) = j • x. 

Plugging this ansatz into (12.61 ) and ordering the powers of e, we find, for the O(e ) 
term: 

9 t 0i + |V0j| 2 = ; <^(0,x)=j-x. 
The solution is given explicitly by 

(2.8) <^(t,x) = j-x-t\j\ 2 . 

The amplitude aj is given by the 0(e l ) and is given by equation (12.1b after project- 
ing the wave along the oscillation e l ^*/ £ according to the the set of resonant phases 
given by Ij (Eqn. (12.2I )). By doing so, we have dropped the oscillations of the form 
e i(k-x-u)t)/e w j t j 1 u | fc| 2 , generated by nonlinear interaction: the phase k ■ x — uit 
does not solve the eikonal equation, and the corresponding term is negligible in the 
limit e — > thanks to a non-stationary phase argument. 

Proposition 12. H is a simple corollary of the following result that is established in 
ll6l . We sketch the proof in Appendix lAl 

Proposition 2.4. Let (a>j) 6 £ 1 (Ij d ), and v e be defined by (12.71) and (12.8I ). Then 
there exists £q{T) > such that for < e ^ £o(T), the exact solution to (12.61) 
satisfies u £ € C([0, T]; W), where T is given by Lemma |23] In addition, v e 
approximates u e up to 0(e): there exists C independent of e e]0, £o(T)] such that 

sup \\n e (t) -v £ (t)\\ w ^ Ce. 

2.3. Link with normal forms. Viewed as an infinite dimensional Hamiltonian 
system, (11.31 ) can also be interpreted as the equation associated with the Hamilton- 
ian 

H e (u,u) = H + eP := ^2 \j\ 2 \uj\ 2 + ^ u k u m U£Uj. 

j&L d k+m=j+e 

that is iuj = dujH(u, u), see for instance the presentations in [llEGl and lfT3l . 

In this setting, the Birkhoff normal form approach consists in searching a trans- 
formation t(u) = u + 0(eu 3 ) close to the identity over bounded set in the Wiener 
algebra W, and such that in the new variable v = r(u), the Hamiltonian K £ (v, v) = 
H £ (u, u) takes the form K £ = Hq + sZ + e 2 R where Z is expected to be as simple 



6 



R. CARLES AND E. FAOU 



as possible and R = 0(u ). Searching r as the time t = e flow of an unknown 
Hamiltonian we are led to solving the honiological equation 

{H ,x} + Z = P, 

where { • , • } is the Poisson bracket of the underlying (complex) Hamiltonian 
structure. Now with unknown Hamiltonians x(u,u) = J2k m j £XkmejUkU m ueUj 
and Z(u, u) = J2k mji Zkm,ejUkU m ueUj, the previous relation can be written 

(\k\ 2 + \m\ 2 - \j\ 2 - \l\ 2 )Xkrnlj + Z kmtj = P km £j, 

where 

f 1 if k + m-j-£ = 0, 

Pkmij S 

[ otherwise. 

The solvability of this equation relies precisely on the resonant relation \k\ 2 + 
\ m \ 2 = \j\ 2 + l^| 2; F° r non resonant indices, we can solve for Xkmij and set 
Zkmij = 0, while for resonant indices, we must take Z km £j = Pkmij- Note that 
here there is no small divisors issues, as the denominator is always an integer (0 or 
greater than 1). 

Hence we see that up to 0(e 2 ) terms as long as the solution remains bounded in 
W, the dynamics in the new variable will be close to the dynamics associated with 
the Hamiltonian 

K[(u,u) = H + eZ := ^ \j\ 2 \uj\ 2 + e- ^ u k u m u e uj. 

j€Z d k+m=j+£ 

|fc| 2 +M 2 =bf+KI 2 

At this point, let us observe that Hq and Z commute: {Hq, Z} = 0, and hence 
the dynamics of K\ is the simple superposition of the dynamics of Hq (the phase 
oscillation ( 12.81 )) to the dynamics of eZ (the resonant system d2.ll )). Hence we 
easily calculate that v e (t, x) defined in ( 12.31 ) is the exact solution of the Hamiltonian 
K\ . In other words, it is the solution of the first resonant normal form of the system 
Oh 

The approximation result can then easily be proved using estimates on the re- 
mainder terms (that can be controlled in the Wiener algebra, see lfT3l ). in com- 
bination with Lemma 12.31 which ensures the stability of the solution of K\ and a 
uniform bound in the Wiener algebra over a time of order 1/e. 

3. An iterative approach 

We now turn to the analysis of the resonant system d2.ll ). The main remark for 
the forthcoming analysis is that new modes can be generated by nonlinear interac- 
tion: we may have aj ^ even though aj = 0. We shall view this phenomenon 
from a dynamical point of view. As a first step, we recall the description of the sets 
of resonant phases, established in (HI in the case d = 2 (the argument remains the 
same for d ^ 2, see O): 
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Lemma 3.1. Let j 6 Z, d . Then, (k, £, m) € i, precisely when the endpoints of 
the vectors k,£,m,j form four corners of a non-degenerate rectangle with I and 
j opposing each other, or when this quadruplet corresponds to one of the two 
following degenerate cases: (k = j, m = I), or (k = £, m = j). 

As a matter of fact, (the second part of) this lemma remains true in the one- 
dimensional case d = 1. A specifity of that case, though, is that the associated 
transport equations show that no mode can actually be created Q. The reason is 
that Lemma I3TT1 implies that when d = 1, (12.11 ) takes the form iaj = Mja,j for 
some (smooth and real-valued) function Mj whose exact value is unimportant: if 
aj(0) = 0, then aj(t) = for all t. In the present paper, on the contrary, we 
examine precisely the appearance of new modes. 

Introduce the set of initial modes: 

J = {j € 7L d | a 3 + 0}. 

In view of (12- lb . modes which appear after one iteration of Lemma I3TT1 are given 
by: 

J 1 = {ieZ rf \Jo | %(0)^0}. 

One may also think of J\ in terms of Picard iteration. Plugging the initial modes 
(from Jo) into the nonlinear Duhamel's term and passing to the limit e — > 0, J\ 
corresponds to the new modes resulting from this manipulation. More generally, 
modes appearing after k iterations exactly are characterized by: 

r j'-i d fe i 

Jk = \j eZ d \\Jj e \ % (0)^0L 

4. A PARTICULAR DYNAMICAL SYSTEM 
We consider the initial datum 
(4.1) u (x) = 1 + 2cos2;i + 2COSX2 = 1 + e ixi + e~ ixi + e iX2 + e~ iX2 . 
The corresponding set of initial modes is given by 

Jo = {(0,0), (1,0), (-1,0), (0,1), (0,-1)} x {0 zd - 2 }. 



It is represented on the following figure: 
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In view of Lemma 13- 1 L the generation of modes affects only the first two coordi- 
nates: the dynamical system that we study is two-dimensional, and we choose to 
drop out the last d — 2 coordinates in the sequel, implicitly equal to Z d-2. After 
one iteration of Lemma [37X1 four points appear: 



,h = {(1,1), (1,-1), (-1,-1), (-1,1)}, 



as plotted below. 



The next two steps are described geometrically: 



As suggested by these illustrations, we can prove by induction: 

Lemma 4.1. Let pel 

• The set of relevant modes after 2p iterations is the square of length 2 P 
whose diagonals are parallel to the axes: 

2p 

:=\JJi = {Uuj 2 )\\jl\ + \j2\^y}. 

1=0 

• The set of relevant modes after 2p + 1 iterations is the square of length 
2 P+1 whose sides are parallel to the axes: 

2p+l 

M (2 P+ i) ._ \J j e = {(jl)j2 ) | maxdjil, |j 2 |) < 2P} . 

1=0 

After an infinite number of iterations, the whole lattice 1? is generated: 

\jM {k) =1? x {rv- 2 }. 
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Among these sets, our interest will focus on extremal modes: for pGN, 

A/? P) := {{h,32) G {(0,±2f),(±2P,0)}}, 

A4 (2P+1) := {(J1J2) G {(±2P,±2P),( T 2 p ,±2^)}}. 

These sets correspond to the edges of the squares obtained successively by iteration 
of Lemma l3.1l on Jo. The set AC defined in Theorem 11.11 corresponds to 

The important property associated to these extremal points is that they are gener- 
ated in a unique fashion: 

Lemma 4.2. Let n ^ 1, and j G A/* . There exists a unique pair (k, m) E 
_/V( n_1 ) x 7\A n-1 ) such that j is generated by the interaction of the modes 0, k and 
m, up to the permutation ofk and m. More precisely, k and m are extremal points 
generated at the previous step: k,m G A/"i™ • 

Note however that points in A/"J n ^ are generated in a non-unique fashion by the 
interaction of modes in Z rf . For instance, (1,1) G J\ is generated after one step 
only by the interaction of (0, 0), (1,0) and (0, 1). On the other hand, we see that 
after two iterations, (1,1) is fed also by the interaction of the other three points in 

M\ '■ , (—1, 1), (—1, —1) and (1, —1). After three iterations, there are even more 
three waves interactions affecting (1,1). 

Remark 4.3. According to the numerical experiment performed in the last section, 
it seems that all modes — and not only the extremal ones — receive some energy 
in the time interval [0, T/e]. However the dynamics for the other modes is much 
more complicated to understand, as non extremal points of AA n+1 ) are in general 
generated by several triplets of points in Af^ . 



5. Proof of Theorem 11.11 

Since the first part of Theorem ll.ll has been established at the end of $2j we now 
focus our attention on the estimates announced in Theorem ll.il 

In view of the geometric analysis of the previous section, we will show that can 
compute the first non-zero term in the Taylor expansion of solution a m (t) of (12.11 ) at 
t = 0, for m G A/"*. Let n ^ 1 and j G N* n \ Note that since we have considered 
initial coefficients which are all equal to one — see (14.11) — and because of the 
symmetry in (12.11) , the coefficients aj(t) do not depend on j G Af* n ^ but only on n. 

Hence we have 

y.a(n) Aa{n) . Wn)+1 fl ja(n.)+l„ . 



aj (t) = , ~ ,7 (0) + — -rr- / (1 - ' (et)de, 

JX ' a{n)\ dt a ( n ) a{n)\ J K ' dH n ) +1 

for some a(n) G N still to be determined. 
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First, Eqn. (12.41 ) in Lemma 1231 ensures that there exists Co > independent of 
j and n such that 

f a(n)+l pi ja(n)+l . 

rj {t) = t —— (i-e)»w l . *M ot)de 

3K ' a(n)\ J Q k ' dH n ) +1 

satisfies 

(5.1) \r 3 (t)\ < (C t) a ^ +1 . 
Next, we write 

(5.2) aj(t) = c(n)t a W +rj(t), 

and we determine c(n) and a{n) thanks to the iterative approach analyzed in the 

previous paragraph. In view of Lemma l4~2l we have 

ihj = 2Xc(n - i)H 2a ( n ~V + o ( i 2Q ("- 1 )+ 1 



where the factor 2 accounts for the fact that the vectors k and m can be exchanged 
in Lemma [42] We infer the relations: 

a(n) = 2a(n - 1) + 1 ; a(0) = 0. 

era) = -2iA — r r ; c(0 = 1. 

v ; 2a(n - 1) + 1 v ' 

We first derive 

a(n) = 2 n - 1. 
We can then compute, c(l) = — 2iA, and for n ^ 1: 

(2A)^=o2 fc . (2A) 2 " +1 - 1 



c(n + 1) = i- 

mti (2 fc - if™ res ^ - iy 

We can then infer the first estimate of Theorem ll.il by Proposition [23J there exists 
C independent of j and e such that for < e ^ £q(T), 

T 

\uj(t) - aj(et)\ ^ Ce, < t < — . 

We notice that since for j £ JV* , |j| = 2 n / 2 , regardless of the parity of n, we 
have a(n) = \j\ 2 — 1. For j £ A/"*, we then use ( 15.21 ) and d5.ll ), and the estimate 
follows, with Cj = c(n). 

To prove the last estimate of Theorem 11.11 we must examine more closely the 
behavior of c{n). Since 2 k - 1 ^ 2 k for k ^ 1, we have the estimate 

22 n+1 -l 

|c(n + 1)l> 2 E£ife2"+i-» - 

Introducing the function 

"+1 x _ x n+l 

/n+i (») = V x k = — x, x^l, 

^— ' 1 — x 

fc=i 



ENERGY CASCADES FOR NLS 



1 1 



we have 



n+l x.x 



n — 3, 



and the (rough) bound 



\c(n + l)\ ^ 2 2 " +1 - 1 - 27l+2 +"+ 3 



-2™ +1 +n+2 



^ 2" 



We can now gather all the estimates together: 

\uj(t)\> c(n)(et) a(n) - (C et) a{n)+1 - Ce 

2™-l 



1 et 



(Coety 



Ce 



(5.3) 



1 /et 



2 n -l 



(1 - (2C ) 2 "et) - Ce. 



To conclude, we simply consider t such that 

2 n -l 



(5.4) 



e 7 , that is t 



£ l-7/a(n) ' 

Hence for the time t given in (15.41 ). since a(n) = \ j\ 2 — 1, we have 

(2C7 ) 2?l et = (2C ) li|2 e^l J l 2 " 1 ) 

= expMj| 2 log(2C ) 
Assuming the spectral localization 

1 



7 



1 



log 



\j\ ^ Q log 



we get for e small enough 



(2C y et < exp a' log - log(2C ) - 4, log - 



20 



1 



1 



1-20 N 



The argument of the exponential goes to — oo as e — » provided that 

7 > and 8 < -, 

in which case we have 1 — (2Co) 2 "et > 3/4 for e sufficiently small. Inequal- 
ity (15 -3b then yields the result, owing to the fact that Ce is negligible compared to 
e 7 when ^ 7 < 1. 

Finally, we note that the choice ( 15.41 ) is consistent with et € [0, T] for e ^ eo, 
for some eo > uniform in j satisfying the above spectral localization, since 

l-20\ 



?7 / a ( n ) 



e Q ( n ) 



7 



< exp ^ log - 



cr 



1 



0. 
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6. Numerical illustration 



We consider the equation (11.11 ) in the defocusing case (A = 1) on the two- 
dimensional torus, d = 2. We take u(0,x) = 5(1 + 2cos(a;i) + 2cos(x2)) with 
5 = 0.0158. With the previous notations, this corresponds to 1/e = 5~ 2 ~ 4.10 3 . 
In FigureQ] we plot the evolution of the logarithms of the Fourier modes log \uj (t) \ 
for j = (0, n), with n = 0, ... ,15. We observe the energy exchanges between the 
modes. Note that all the modes (and not only the extremal modes in the set A/"*) 
gain some energy, but that after some time there is a stabilization effect (all the 
modes are turned on) and the energy exchanges are less significant. 



Figure 1 . Evolution of the Fourier modes of the resonant solu- 
tion in logarithmic scale: uq(x) = 1 + 2cos(xi) + 2cos(x2) 

In contrast, we plot in Figure [2] the solution corresponding to the initial value 
u(0,x) = <5(2cos(xi) + 2cos(x2)) with the same 5. In this situation, no energy 
exchanges are observed after a relatively long time. Note that in this case, the initial 
data is made of the 4 modes {j £ 1? ||j| = 1 } forming a square in 1? . This set 
is closed for the resonance relation, so no energy exchange is expected in the time 
scale 0(1/ e). We notice that the solution of (12.11 ) is given explicitly by aj(t) = 
f° r bl 7^ 1' an d a j(t) = exp(9it) for \j\ = 1. 

The numerical scheme is a splitting time integrator based on the decomposi- 
tion between the Laplace operator and the nonlinearity in combination with a 
Fourier pseudo-spectral collocation method (see for instance fi~8l and jlTl Chap 
IV] for convergence results in the case of dl.ll T). While the Laplace operator part 
idfU = —Au can be integrated exactly in Fourier, the solution of the nonlinear 
part idtu = \u\ 2 u starting in v(x) is given explicitly by the formula u(t,x) = 
exp(— it\v(x)\ 2 )v(x). The fast Fourier transform algorithm allows an easy imple- 
mentation of the algorithm. The stepsize used is r = 0.001 and a 128 x 128 grid 
is used. 

Note that using the framework of lfl2l [Til , we can prove that the numerical 
solution can be interpreted as the exact solution of a modified Hamiltonian of the 
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Figure 2. Evolution of the Fourier modes of the nonresonant 
solution in logarithmic scale: uq{x) = 2cos(xi) + 2cos(x2) 



form 

jeB K (k,m,j,e)eB K PV fcm ^ y 

k+m-j-l£KI? 

where oo k mij = |/c| 2 + |m| 2 — 1^| 2 — |m| 2 and Bk the grid of frequencies j = (ji,j2) 
such that ji and j'2 are less than K/2 = 64. Note that this energy is well defined as 
TOJkmtj is never a multiple of 2tt, and that the frequencies of the linear operator of 
this modified energy carry on the same resonance relations (at least for relatively 
low modes). This partly explains why the cascade effect due to the resonant system 
should be correctly reproduced by the numerical simulations. 

Appendix A. Sketch of the proof of Proposition I2.4I 
By construction, the approximate solution v e solves 

ied t v £ + e 2 Av e = Ae|v £ |V + Xer £ , 

where the source term r £ correspond to non-resonant interaction terms which have 
been discarded: 

j£Z d {k,i,m)^Ij 

We write 

<f) k (t, X) - 4>l(t, x) + 4> m (t, X) = Kk,£,m • X - UJk,£,mt, 

with K kt i jTn € Z d , u kitjm e Z and \n k ,e, m \ 2 ¥= u kj(>m , hence 
(A.l) ||Kfc^ jm | 2 — u} k ,e, m \ ^ 1- 

The error term w e = u £ — v £ solves 
ied t w £ + e 2 Aw £ = Ae (|w e + v £ | 2 (w £ + v e ) - |v £ |V) - \er £ ; wf f=0 = 0. 
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By Duhamel's principle, this can be recasted as 

w £ (t) = -i\ f* e i£{t ~ s)A (|w £ + v £ | 2 (w e + v e ) - |v e | V) (s)ds 
Jo 

+ iX f e i£{t - s)A r £ {s)ds. 
Jo 



Denote 

rt 



R*(t) = [ e i£( - t -^ A r £ (s)ds. 
Jo 



Since W is an algebra, and the norm in W controls the L°°-norm, it suffices to 
prove 

ll-R e ||L°°([o,T];WO = 0(e). 

We compute 
where 

, ,, s /"* , x_ , v , s f .Kk,£,rn ■ % + \Kk,e,m\ 2 S — OJk,e,mS\ , 

Ok,i,m[t,x) = J ak{s)ae(s)a m (s) exp U ■ — Ids. 

Proposition I2.4l then follows from one integration by parts (integrate the exponen- 
tial), along with (IA. 1 b and Lemma 1231 

Acknowledgements. The authors wish to thank Benoit Grebert for stimulating dis- 
cussions on this work. 
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